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Abstract 

This is an extended version of the talk given by the Author at the 
40th Symposium on Mathematical Physics held in Torun, Poland, June 
25-28, 2008. We review the methods of canonical quantization of free par- 
ticle motion on curved submanifolds considered as a system with second 
class constraints. The work is based on our previous articles, [I] and [2]. 
However, some new results are also presented. 

1 Introduction 

We consider the problem of quantum motion in curved spaces. It is well- 
known that in the case of Euclidean spaces the correct quantum Hamiltonian is 
H = — ^-A. Podolsky [3] in 1928 proposed that for arbitrary spaces it should 

be replaced by H = — ^Alb with Alb being the Laplacc-Bcltrami operator. 
This postulate is a direct and geometrically clear generalization of the dynamics 
in Euclidean spaces. But if one wants to get the theory by some canonical pro- 
cedure, he encounters a severe problem. For any given classical theory there is 
an infinite number of quantum theories with a proper h — > limit. Quantization 
is not unique. In Euclidean spaces Dirac recipe in Cartesian coordinates yields 
experimentally correct result for the theories which we usually have in theoret- 
ical physics. However, in curved spaces we do not have a notion of Cartesian 
coordinates and can't make a choice of the theory in this way. Of course, if 
the theory posesses a large enough symmetry it can be sometimes completely 
defined by the symmetry requirements directly at the quantum level, see g]. 
Nevertheless, we think it is very instructive to study the properties of classical 
quantization methods for constrained systems since theories with constraints 
are so important in modern physics. In particular, a possible solution for our 
problem is to embed the space under consideration into some Euclidean space 
and to quantize the new theory as a theory with second class constraints. As 
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we shall see below, the results depend both on the choice of embedding and on 
the method of quantization. 

In Section 2 we describe the Dirac approach in the case of codimension 

1 surfaces first, give an explicit operator realization of it and show that the 
Dirac procedure is ambiguous. The quantum Hamiltonian depends even on a 
particular form in which the equation of surface is presented; a natural geometric 
way of fixing this freedom is explained. We also compare these results [1] with 
those obtained by general relativity inspired methods in [5] and discuss the 
situation in higher codimcnsions. 

In Section 3 we present the method of converting the second class constraints 
into the first class ones by adding some new degrees of freedom. An error 
contained in [T] is corrected. 

In Section 4 we review the thin layer quantization. In this approach a particle 
moves between two equidistant infinite potential walls [6] or it is subject to some 
potential force which in a proper limit makes it moving strictly along the surface 
[7J. In higher codimensions the method becomes quite involved and technical, 
but it also exhibits interesting features [8] related to a remarkable appearance of 
nontrivial gauge structures (for quantum motion in submanifolds with non-flat 
normal bundles) first reported in [351 EH] • 

2 Dirac quantization 

Suppose we have a theory with 2N constraints (f> a , a — 1,2,..., 27V. These 
constraints are said to be of the second class if det{</> a , 4>b} ^ even in a weak 
sense [5] in Dirac terminology (on the constrained surface (p a — 0, Va). In 
particular, it means that, unlike the first class constraints, they do not form 
a closed algebra with respect to the Poisson brackets. In such a situation the 
standard replacement of Poisson brackets by commutators doesn't work as it 
would contradict {<pa,4>b} ^ inequality, and more complicated procedures are 
needed. One possible way out is to introduce the Dirac brackets: 

2N 2N 

if, g}v = {/, 3} - E E{/< <MA„ 6 {&, g}, (i) 

a=l 6=1 

where A a {, is the matrix inverse of {4> a , 4>b}- Now {4>\, §i}d — and we can 
introduce the commutators in a usual way. Dirac bracket is degenerate and does 
not define any symplectic manifold but it can be regarded as a Poisson structure 
obtained by factorization of original Poisson bracket algebra over motions in 
unphysical direction [T], see the expression ([7]) for the momenta operators below. 
As a simplest example, one can consider a free particle motion on (n — 1)- 

n 

dimensional sphere, ^2 xf — R 2 , in n-dimensional Euclidean space. It can be 
considered as a system with two second class ({</>i, fa} = 2 x ^ 0) constraints 
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n 

1 =£a;?-JZ a = O J (2) 

i=l 

n 
i=l 

where pi are canonical momenta. We impose the cf>x condition as a primary 
constraint directly at the Hamiltonian level. In the Lagrange formulation we 
also have a Lagrange multiplyer which would generate one more pair of second 
class constraints (vanishing of its canonical momentum as a primary constraint 
and some consistency condition involving the multiplyer itself). For the physical 
phase space it influences neither the Poisson structure nor the physical Hamil- 
tonian [11) . Thus we prefer not to increase the number of unphysical variables 
more than needed and work purely within the Hamiltonian mechanics. (Note 
that every new constraint enters the Dirac generalized Hamiltonian multiplied 
by its own new arbitrary function of time [S].) 

A natural question to be asked at this point is whether one could quan- 
tize the theory without any classical rearrangement before and then impose the 
constraint 4>\ at the quantum level. It would be problematic because the sec- 
ondary constraint <f>2 (i-e. consistency condition concerning the time derivative 
of 4>i ) shows up which has a non- vanishing Poisson bracket with <j>\ . The prob- 
lem would be to find a natural selfadjoint restriction of the Hamiltonian to the 
physical space. It is possible in principle but complicated and requires quite 
some accuracy. We refer the reader to the Ref. [T2] and proceed with the Dirac 
method. 

A simple calculation according to |T]) shows that [TU] 

{xi,Xj}T> = 0, (4) 

{xi,Pj}v = Sij ~ -z^T> ( 5 ) 

x 

{Pi,Pj}l> = —{Pi x 3 ~ PjXi)- (6) 
X 

This algebra can be satisfied by very simple (usual) coordinate operators x% — 
Xil and standard differential operators of momenta from which the normal dif- 
ferentiation is subtracted pQ: 

-iHv^-ih(v-^(^-v))=f. (7) 

The problem is that pi are not selfadjoint. But at the sacrifice of Leibnitz rule 
we can introduce new selfadjoint momenta: 



2 1/* 4> A .71—1 Xi f 

Pi = ~iPi + Pi) = Pz + ih^r- ■ L 
z z x 



3 



The primary constraint @ defines the space of physical states and the secondary 

n 

one ([3]) acquires the form of identity <j>2 — X (%iPi + i&iPi) ) = 0- The resulting 
Hamiltonian pQ 

~2^ Pl 2 LB 8i? 2 

contains a typical quantum potential Vq D ^ — - ■ We should stress that 

this result can be obtained in a purely algebraic manner without any kind of 
explicit operator realization [10 . 

Note also that the same procedure may lead [1] to Podolsky theory if one 

n 

takes our initial definition ([7]) for pi and Hamiltonian — i X PiPi which 

i=l 

equals — ^-Alb for the physical sector functions. The quantum potential is 

(v) 

zero: Vj r> = 0. Thus one preserves an important property of momenta op- 
erators, the Leibnitz rule, so that they are differentiations on the algebra of 
smooth functions. These operators are not selfadjoint and can't represent ob- 
servables. But in any case they do not have any clear physical meaning being 
projections of generators of motions along the coordinate lines of n-dimcnsional 
flat space, which are somewhat esoteric for an observer living on the sphere. 
Natural observables on the sphere are generators of SO(n) rotations, and they 
are selfadjoint (proportional to i[pi,pj]). 

2.1 Arbitrary codimension 1 surfaces 

The free motion on a codimension 1 surface f(x) — can be obtained [TT] by 
an obvious modification of © and ([3]): 

0! = f(x) = 0, (8) 

n 

fa EE Yjfrf)* = 0. (9) 



If 



V/ 



7^ at the physical surface (as we assume throughout the paper), 

these constraints are of the second class because {</>i,</>2} = fv/) ■ Again we 
introduce the Dirac brackets by the prescription (fTJ) and get 

{xi,Xj}v = 0, (10) 

{XhPjjv = kj — 2 , (11) 

V/ 



n 

{p^P 3 }v = -j— T J2 ((^/)(^-/) - (dif)(dy)) Pk . (12) 

V/) k=l 
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One can use the following operators [T] for the quantum description: xi = Xil, 



Pi 





d 


V/ 


dXj 



as non-selfadjoint momenta and 

ih ■ 



( 



2 f-f cb , 

— 1 \ J 



V 



v/ 



(13) 



for the selfadjoint counterparts. (Here we implement the factorization over 
unphysical motions again.) The operator ordering problem is relevant only 
for the momenta commutators, and it is solved by our explicit choice of the 
operators as follows: 



\pi,Pj 



V/ fc=1 



\pM = VE 



fe=l 



— — Pk + Pk 



V 



(v/)' 



V/ 



/ 



We 



n n / ^ ^ \ 

have the identity Y^(^if)Pi = or £ [(@if)Pi +Pt(^*/)J = for the 
i=i i=i ^ " ' ' 

secondary constraint. And the physical sector is defined by the primary one: 

Vphy. = tp(x)S(f(x)). 

For non-selfadjoint momenta the Hamiltonian reads 



H {V) = \±P\P, = -^U- 



d 

dn 



div(~u 



d 
dlt 



where A is the Laplace operator in the Euclidean space and n 



_ v/ 



-±r-T is a 

V/ 



unit vector normal to the surface (JSj> . In the selfadjoint case the Hamiltonian 
= h j>2 Pi = + Vq V \x) contains also a quantum potential 



^ = -yE E 



8=1 \ j = l 



dxj 



E 



9xi 



V 



E 

fe=i 



(v/) 

(gi/Xgfc/) a 
dxk 



(V/)' 



E 

V =1 



V/ 



(14) 



5 



Some details of rather straightforward calculations can be found in pQ. Note 

that in terms of the normal vectors it — tSA- one can easily write the potential 

|v/| 

(Q3|) down in the following form: 



If 
4 



Vq V) = x ( \ ( diUi ) + H dl ( n k d kni) + nin k n m dl m ni ) (15) 

\ \ i / i,k i,k : m 



using the obvious relations Y fiidkTii = and Y n i°\m n i = ~ S (°m n i) (^fc^i)- 

i i i 

Unfortunately both Hamiltonians, and H^ v \ are ambiguous; they take 
different values for those functions which represent one and the same surface. 
(And the problem exists even for spheres.) Indeed, any surface can be repre- 

n-l 

sented by its tangent paraboloid at some point: f(y) = y n — \ kaVa + ®{Va) 

a=l 

where y a are Cartesian coordinates. (It is not a priori obvious that this accu- 
racy is enough for calculating the quantum potential but one can easily check 
that in this case it is, see [T].) Then (fT4")) gives 

n-l \ 

in the vicinity of the point y = 0. For a sphere the principal curvatures are 
k a = ji and at the chosen point we have V q — g^a-^ which differs from 
our previous result (and that of |10j). So, the Dirac recipe is ambiguous. To 
fix the freedom, we propose the following (geometrically natural) choice of the 
function f(x): up to the sign it should be equal to (some function of) the 




distance from the surface / = 0. After that we have 



V/ 



= 1 and = dkn t , 



^rifc(?fc7ii = where n^ = dkf- Then a simple calculation [1] shows that the 

k 



quantum potential on the f{x) = surface is V q = j E • For spheres 



a=l 

it yields the previous result. The kinetic part of the Hamiltonians for our choice 
of unit normals 

A - ( -JL ) - div(^) • Tpr — £^lb 
\on J an 

equals to the Laplace-Beltrami operator on the physical surface pQ. But in the 

general case (when the unit normal vector ~n = \ would not be orthogonal 

to the surfaces f(x) = const ^ 0) this result would not be true. The vector 
n will change its direction while moving apart from the initial surface, and the 
second normal derivative (^=?-) would add some extra (first order differential) 
term to Alb- 

Let's consider a simple illustration. For a circle in a plane we would use 
f(x,y) = \y\ — Vl — x 2 (instead of f(x,y) = x 2 + y 2 — 1) and approximate it 
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near the (x,y) = (0, —1) point by a parabola f(x) = y — ^- + 1 = 0. We have 
n x = ~—==2 and n y = T/fTjf ■ The selfadjoint momenta can be easily found 

as 

ih f d d \ ihx 

= — 9 \ + x - 



1 + x 2 \dx dyj (1 + x 2 ) 
ihx f d d\ ih (x 2 - 1) 



l + x 2 \dx dyj 2(1 + a; 2 ) 2 
It leads to the Hamiltonian 

_ h2 ( 1 ( d , „ 9 \ 2 2x ( d , _ d \ \ , h2 ( 3 - 2x2 - 5x4 



2yi + x 2 \dx dyj (l + x 2 ) 2 \dx dyj J 8{l + x 2 ) 4 

iV q — ^g- for the quantum potential at 
introducing the tangent derivative S=r = , 1 ■ I - 



which gives V q = for the quantum potential at x — > (instead of ^-). And 




2.r 



LB 



(i + .x 2 ) 3/2 at / 

for the kinetic energy operator with Alb = 

d t 

2.2 Relation to GR-like methods 

Some time ago this problem was tackled in Ref. 5J by methods typical to 
general relativity. The quantization was performed in Cartesian coordinates but 
a curvilinear coordinate system was also used. One of the coordinates qo was 
chosen to be the value of the function / and the others had to be orthogonal to it. 
Then the authors of [5] defined the curvilinear momenta p M in terms of Cartesian 
ones. The Cartesian commutators (obtained from Dirac brackets) implied the 
commutation relations for the curvilinear operators [q u ,p^\ — ih (<5^ — n^n^ 
and set po — 0. After doing this one can arrive at [5] 

n — 1 

where gij is the metric on the physical surface and the quantum potential is 
equal to our result for the special choice of the equation of surface, ^ (divT?) 2 . 

The authors interpreted the kinetic part of the Hamiltonian as the Laplace- 
Beltrami operator on the hypersurface / = 0. It would indeed be true if the 
momenta were the standard ones for the hypersurface, but they were defined 
to be symmetric in the ambient space which is not the same. Up to the 
setting the normal differentiation to zero, they would be the standard sym- 
metric momenta operators for the curvilinear coordinate system in the am- 
bient Euclidean space. It amounts to the difference between -r/idig 1 ^ 4 and 
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^7^9; ! ■ If the function / depends only on the distance from the surface 



/lv/1 

(di\ V f\ — 0) then it makes nothing and our results coincide. What changes 
if we have another function /? A linear differential operator ^- Y2 ^Tjp S y &j 

gets added to — ^-A^b and two additional terms appear in the quantum po- 
tential: — ^ • — h ipdiv V [^/(^ ■ AH the operations are related to the 
hypersurface / = 0. In particular, if the gradients are taken in some Cartesian 
coordinates of the ambient space, the differentiations should be projected to 

the hypersurface. In this coordinates the vector ^lyjp in the tangent space 
of / = would have the following components: j^jj ^di — riiY2 nkdk^j | V f\ = 
Y] J^trdj . And the new terms in the quantum potential can be transformed 

3 

to \ (rijOjUi) — ^g-J2 [ ^rijOjUi in a complete accordance with (TT5|) . 

i,j i \ 3 j 

The differential operator converts to Y2 rik{dkni)di in Cartesian coordinates, 

i.k 

which is exactly what would come out of the second normal derivative (g=s-) 2 
in our approach. 

2.3 Dirac quantization in higher codimensions 

Going to higher codimensions (the surface = for a = 1,2... N) compli- 
cates the things considerably because the Poisson brackets of Y^Pi^if and 

i 

Y^Pi^if (the secondary constraints) do not vanish even if y/^ ■ y = 0. 

i 

It affects the momenta commutators. Nevertheless, if the normal bundle of the 
submanifold is flat and we can (locally) choose the set of functions / satisfying 

the conditions y/ (a) = 1 for all a and y/ (a) ' V/ (fc) = for a ^ b then 
quite bulky but straightforward calculations show that the generalization of the 
previous results is very simple: 

^ N n „ \ n Q / N \ 



h2 N ta 2 



a=l 



It gives the extrinsic mean curvature squared for the quantum potential. 

And in the general case of higher codimensions, the straightforward approach 
equations become almost untractable. But the structure of terms still remains 
quite understandable. For example, in the case of codimension 2 one can check 
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that the momentum operators can be obtained by projecting the vectors to 
the subspace orthogonal to both v/^ and V according to 




At this point the methods of previous subsection [5] become very useful. In 
[I"3] it has been shown that these methods lead to formally the same result 
as for the codimension 1 case (Laplace-Beltrami operator and extrinsic mean 
curvature squared). As it was above, the interpretation of the kinetic energy 
operator is correct if the normal components of the metric do not influence the 
momenta operators which means that det (\/f^ ' V/^J does not depend on 
tangential coordinates. 



3 Abelian conversion method 

Sometimes in the quantum field theory first class constraints may fail to form 
a closed algebra at the quantum level. Formally, they can be said to become 
second class upon quantization, but it is very bad for the quantum field theory 
with local symmetries as it corresponds to a gauge symmetry breaking. For 
example, this kind of anomaly occurs for the Gauss law in the model of Weyl 
fermions interacting with a Yang-Mills field. In [14] it was proposed to introduce 
some new degree of freedom in the model in order to have new constraints with 
the Abelian algebra and get rid of the anomaly. It resembles the appearance 
of the conformal factor as a new degree of freedom at the quantum level for 
noncritical strings [15] . (Very similar situation occurs in the theory of relativistic 
branes [16] , but it is not so clear what to make out of it in this context and also 
it is not known if the critical dimension exists, see [T71 HU [19] for diverse views 
on the subject.) 

This idea can also be used in quantum mechanics with second class con- 
straints [10] . The ( "Abelian conversion" ) method consists of introducing new 
canonical pair of variables Q, K and first class constraints eri, 02 : {ai, 02} = 
(identically, so that the constraint algebra is Abelian) assuming <j\ = (f>i , 02 = 
<p2 if Q = and K = 0. In our case it would be o\ — f(x) + K and 

o"2 = n ■ p + 



V/ 



V/ 



Q. (Note an error at this point in [JJ where the 
factor has disappeared from 02 •) The next step is to find a new Hamiltonian 
such that H s = H if Q = and K = and {H s ,ai} = {H Sl <7 2 } = 0. The 
physical sector is obtained by setting o\ = er 2 = 0. For free motion on spheres 
this method gave zero quantum potential, see [TUJ [TTJ . Actually, the authors of 



[TOl [TT] had the result of the form Hs — Hg f ai, 02, X) ( x iPk — XkPi) 2 ) which 

i<k 
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due to relation 



i<k 



could be transformed to H, 




phys 



because <rf 



(Ew) 2 if Q = o. 



It was concluded that the Abelian conversion method is preferable because 
it involves no extrinsic geometry in the results. For this and other reasons it 
was used in the beautiful projection operator approach to path integral quan- 
tization of constrained systems [20] aiming at quantizing the gravity |21) . But 
is it possible to generalize the above result to other surfaces? In principle, 
methods which introduce new variables are very strong, see for example [22| 
and references therein. Moreover, one can proceed with a general philosophy 
of Abelian conversion [11) without any need of introducing such weird concepts 
like ghost operators [22]- The question is whether it is possible to do this using 
the very simple geometric form of the new constraints introduced above and 
obtaining some relatively simple and physically sensible results. Let us search 
for Hs = Hs(crx,cr2,x,p) in the previous form 



Hs = H s \ ai,a 2 ,g(x) 




(16) 



designed for getting the pure Laplace-Bcltrami solution. The following relations 
show up: 



{H s ,<r 1 } = -'£ i i 



' dpi 



0, (17) 



{H s ,(r 2 } = 

i 

From (fl8l) we have 



Eft(^ - E 1 v/|) = <>. , is) 



dH s 

oxi t-r* op, 



i.k 



dpi 



y^.PiPk \ ^2 n ]d 3 g(x)(S ik - mn k ) - 2g{x){din k ) 

~ 2 E { 9l V/ ) Q ^Pi ~ n iYl n i P 3 ' j 



If we admit the condition 



V/ 



= 1 (compare with the Dirac method above) 

then the last term disappears and it has a non-zero solution for spheres because 
diUk ~ Sik — n-iUk- But this is not true for arbitrary surfaces. Hence the result 
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of [TU] can't be generalized directly. Nevertheless, using a simple ansatz 



\ i,k i 

with Cik — Cki one can show [1] that in general it is possible to get a quadratic 

in momenta Hamiltonian (not equal to the Laplace-Beltrami operator) by this 

method if we admit the above definition of the function f(x). 

If another function is used then (fl8|) can't be valid identically because the 

last term contains Q. But we can afford having it only in a weak sense. Then 

from cr 2 we determine Q = — il^ p i , and (fTgl) with the ansatz (Til)]) converts into 

v/ 



^ rijdjg(x) (S ik - n;n fe ) + 2g(x) [ m ^ njdjUk - d l n k 



which is probably not as hopeless as it was erroneously stated in [T] due to the 
aforementioned mistake but still can give no guarantee for the existence of a non- 
trivial solution. Another problem of the method is that in this setting there is 
no clear reason for insisting on the Abelian algebra of the first class constraints. 
And if we go to higher codimensions we would really need to modify the method 
somehow, at least by finding a more clever choice of the new constraints because, 
as we already mentioned, the Poisson brackets of ^Pinf^ and ^Pinf^ do not 

i i 

■ i -r — >( a ) — K b ) n 

vamsn even it n ■ n = I). 



4 Thin layer quantization method 

As it was discussed in the Section 2, imposing the second class constraints 
directly at the quantum level is problematic. Nevertheless, we can use a more 
delicate procedure. We can approximate the constrained system by a motion in 
a thin layer around it. In quantum mechanics this approach appeared in [51(7], 
for a deeper discussion see [8] [23] and [24] at the mathematical level. It can 
also be used for classical systems [M] and gives rise to extra potential if we take 
the initial conditions involving normal motions in the thin tube, but coincides 
with the intrinsic description if the initial velocities are tangential. However, in 
quantum mechanics this approach always gives a geometric potential because in 
the quantum realm it is impossible to eliminate the normal motion completely. 

We consider (n— l)-dimensional smooth surface in M n and two infinite poten- 
tial walls at the distance 8 — > from the surface. Free quantum particle moves 
in the thin layer of width 26 between these potential walls. We introduce a 
curvilinear coordinate system in which \x n \ equals the distance from the surface 
to the given point (thus playing the same role as the function / in the refined 
approach to Dirac quantization), and the coordinate lines of x±, . . . , x n -\ are or- 
thogonal to that of x n . We have the boundary condition $>\ x _ s = ty\ x __ 5 = 
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and Hamiltonian H = — 4j-A with A being the Laplace operator, 

n n 

A = E E ~9~ 1/2 dfg 1/2 ~g lk d k = d n 2 + (r 1/2 a„5 1/2 ) d n + A LB , 

i=l k=l 

g a b o 
o 1 

where Alb is the Laplace-Beltrami operator on the surface x n — const. 

The simplest way [5] to obtain the thin layer limit is to consider the tangent 

n-l 

paraboloid of the surface y n = | ^2 ^aVa + ^iUa)* where k a are the principal 

a=l 

curvatures. The unit normal is n a = , kaV " + 0(y 2 ) = k a y a + 0(y^), 



n n = -1 + 0{y 2 a ) and 

n-l 

divlt = J2^ + 0(y a ). (19) 

a=l 

A nearby surface i„ = e can be obtained by taking ~y — ► ~y' = ~y + e~n and 

"U\l+0(y' a 2 ))dy' a n-l 

dy' a = dy a (1 + ek a + 0{y a )). It yields % = = fl (1 + e£: a ) + 

n (i+0(y^))^a <*=l 

a = l 

0(y a ) near the point y = 0. At the line y a = Va = 1, . . . , n — 1 one has 

§ = 1 +^E^+r 2 (fE^ -E^+° (e3) - (20) 

a— 1 y \a— 1 / a— 1 / 

Clearly, the relation (|2"D|) is valid for every point of the surface with its own 
principal curvatures. Following [7] we introduce a new wave function 

X (x) = 

Physically it amounts to 

6 

2 



dV\^{x)Y = J dx n l dS\ X (x)\\ 
\x n \<5 -S 

so that the function J dx n \x(x)\ 2 defines the probability density of finding the 
particle at a given point on the surface. For the lowest energy solution the 
normal motion gives only the factor of cos ^f- , and we easily get [2] 



A*(s) = A LBX (x) + dlx{x) + ( i k 2 a - \ [J2 ka) ) X(x) + 0(x n ). 



a=l 
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At this energy level x( x i, ■ ■ ■ ; x n) — f{xi, ■ ■ ■ , x n -i) cos^f^. After taking 5^0 
limit and subtracting an infinite (proportional to l/<5 2 ) energy we obtain the 
Hamiltonian 

H = ~^ + Y\p\ ka ) ~ 2 2 fc ') (21) 

which contains the quantum potential 

^(&) 2 - 2 !'4 

For 2-dimensional surfaces in M 3 the result of da Costa [7], V q = — ^-(fci — fe) 2 , 

is reproduced; for spheres k a — i and the potential is V 9 = - ^|ji"~ 3 ^ ■ If we 
would use a layer of varying width, some additional effective forces will appear 



4.1 Some remarks and variations 

We could use an appropriate confining potential instead of infinite walls. It 
would lead to the lowest energy level function of the potential V con f(^f) instead 
of cos ^f- and to another infinite energy. Note that we can also embed one 
curved space into another curved space [5]. Moreover, this approach can be 
used for the quantum graphs theory [35] . And it is currently used to describe a 
motion of electrons in nanostructures [SSj [STJ and for the physics of molecules 
[HI HE] • Even before the general consideration of the problem appeared in |6J 
and [7J, some elements of the thin layer approach were successfully used in 
the theory of chemical reactions [29 . Recently, a considerable progress in the 
method has been achieved [41] allowing one to treat a very general type of 
quantum constrained motion, even with relatively large kinetic energies, with a 
full mathematical rigour. 

We should mention that there is one more method of quantization proposed 
by Prokhorov [30j . The motion of a particle is considered as a system with 
two second class constraints but only one condition is imposed on the physical 
sector: P n ^ phys (x) — with P n = —ih-rjjj-^-g 1 / 4 . It means that 

dn (VW^^^ = ^ (22) 

Having solved some task by this method, one should put x n = in the results 
after all the differentiations over x n are performed. Due to l[22|) the probability 
to find a particle at the distance \x n \ from the surface does not depend on the 
value of x n , and we choose one value we need. (For Prokhorov's view see [50].) 
This method gives the same results [5] as the thin layer approach due to a very 
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simple physical reason. The lowest energy level wave functions (in the model 
with two infinite potential walls) have nodes at x n — ±5 and the bunch at 
x n = 0: dnX = or, equivalently, P n ^ = 0. 

One also could be tempted to use the Hamiltonian in curvilinear coordi- 

fn-l 2 \ 

nates H = | ( Yl T^fey + Pn) w ^ n a simple recipe pi — > pi = —ihg^^-^dig 1 ^ 

followed by the thin layer method. Then the result [31] for ip(xi, . . . , x n ) = 
f(xi, . . .,X n -i) ■ cos ^§f- 

/ . cos ^ _ 1 ^ p| /. cos ^ 1 h 2f 7TX n 

[7E *hhi{x) [7F 2 [7F n 25 ' 

V dS 1-1 V dS V dS 

can be considered as zero quantum potential. However, one should remember 
that quantization in curvilinear coordinates is dangerous because its results usu- 
ally depend on the choice of coordinate system. And the curvilinear momenta 
operators are only symmetric but not in general selfadjoint. And what is more 
important, the operator ordering problem in fe p ^ terms is not solved. It is 
not difficult to deduce the correct ordering for the zero potential theory, but 
this particular ordering is not so natural a priori and can involve quite bulky 
expressions [2]. (See also [32].) 

4.2 Higher codimensions and gauge structures 

In general we can represent a smooth m-dimensional surface in IR™ by its tangent 
paraboloid at a chosen point: 

m m 

S« = 2EE^» + «, (23) 

a=l 6=1 

a = m + 1, . . . ,n with some "curvature coefficients" — . But it turns 
out that this approach is convenient only for quantization on curves (and for 
some other relatively trivial cases like flat 2-torus in R 4 ). 

4.2.1 Quantization on curves 

For 1-dimensional manifolds (curves) a suitable rotation in the space of y a casts 
([23]) to the form y 2 = \ky\ + 0(yf); 2/3,..., y„ = 0(y\). The unit normal 

vectors are nf } = k Vl + 0{yf), n ( 2 2) = -1 + 0(yj ), nf ] = ... = = 0(yf); 

n 2 (a) = -5 ia + 0(yf) for a > 3. We have ~n [a] li {P] = 6 af3 +0(yf), and after the 
/ n ( ) 

transformation ~y — > ~y = ~]f + £a~n one gets dy[ = (1 + e 2 k + 0(yi))dyi 



a=2 

and dy' a = (1 + 0(y\))dy a for a > 3. We introduce a new curvilinear coordinate 
system near the curve in which x\ is just the length along the curve and the 
hypersurfaces of constant x\ are the cross sections of its tubular neighbourhood. 
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And in a given cross section any point ~~r has other n — 1 coordinates defined 
by x a = n ■ r . In this coordinate system 



Hk = 



{l + x 2 kf 
I 



and A = A c + A„ + ( i+^^l + x 2 k)) d 2 = A c + A„ + jj^d 2 where 



n 

A c is Laplace-Beltrami operator on a curve x a = const and A„ = d 2 , 
is Laplace operator in a hyperplane x\ = const. Then for a wave function 



= ^/T+H^k ^>(x) in a thin layer ^2 x^ ^ S 2 we obtain [2] 

a=2 

A*(x) = A— ^=L= = A cX {x) + A nX (x) + ^-x(x) + 0(x a ). 
^/l + x 2 k 4 

After subtracting an infinite energy due to A n \(x) it yields the quantum po- 
tential V q = — \k 2 as in [?]■ At this point we can also see that the higher 
codimensional thin layer problem is not reducible to a step-by-step decreasing 
of the physical space dimension. Indeed, for a straight line in R 3 we obviously 
have zero quantum potential. But if we first consider a cylinder of radius R and 
then restrict it further to the line, the quantum potential would be — g^r- 

There is one subtlety in the above discussion. If the curve has a torsion, our 
coordinate system rotates around it. If one attempts at describing this motion 
in a non-rotating coordinate system he gets new terms in the Hamiltonian which 
correspond to rotations around the curve [33 . Locally these descriptions are 
equivalent, but globally for a closed curve the rotating coordinate system may 
not exist and one can get global phases out of it [33] . (Global effects may also be 
relevant for nontrivial normal bundles, for example for a motion on the Mobius 
strip, [24] .) Note that in our consideration we need either to use a thin tube with 
spherical cross sections or to rotate the cross section around the curve together 
with its Frenet frame. For more general discussion of twisting the confining 
potential see [8]. 



4.2.2 Higher dimensions and geometry of normal bundles 

In general we can try to use a similar construction; but such coordinate systems 
which eliminate all the rotations from the Hamiltonian do not exist for subman- 
ifolds with non-flat normal bundles. It was first noticed by da Costa [33] who 
proved that it is not always possible to find a smooth family of normal vectors 
with the properties required for separation of normal and tangential motions. 
Geometrically we need to understand the structure of a tubular neighbourhood 
of the physical submanifold which is naturally related to a small portion of the 
normal bundle corresponding to its embedding into the ambient space. Then 
one has to define the standard mathematical notion of the normal connection 
in the normal bundle [U [35] . We will not discuss it here in any detail, but basi- 
cally it amounts to taking the normal projections of the ambient space covariant 
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derivatives. From the ambient (Euclidean) space viewpoint parallel transports 
according to this connection involve rotations including those which go around 
the physical subspace. For example, in the tubular neighbourhood of a curve 
the normal connection rotates the normal vectors not only with the normal hy- 
perplanes but also around the curve together with its Frenet frame. Of course, 
on a curve we can exclude this effect by a suitable choice of rotating coordinates 
[2] as in the previous subsubsection. But for higher dimensional (and codimen- 
sional) manifolds the normal bundle can be non-flat and such exclusion would 
be impossible. The simplest examples are the configuration space of the double 
pendulum [33] and the helical surface from [33]. 

Analytically it means that the Laplace operator necessarily contains terms 
with mixed normal and tangential derivatives, d\ a . They vanish at the physical 
surface linearly with the distance from it, bit it is obviously not enough for the 
thin layer approach (unlike the Prokhorov method). Da Costa concluded [34] 
that the thin layer quantization would not work well in this situation. This state- 
ment was repeated in [5] . Strictly speaking, this attitude is not right because all 
the dangerous terms sum up to the angular momenta operators corresponding 
to rotations in normal sections of the thin tube [8] [23] [39l [40] . In the case of the 
lowest energy solutions for a thin tube with a spherical cross section it will not 
influence the resulting theory at all. But if the chosen normal energy level is de- 
generate then a gauge structure will show up [HI [231 EH1 (33 [351 SO] • The simplest 
example (although not very natural from the thin layer quantization viewpoint) 
is a higher energy level of normal motion in the thin tube [23] . In any case, 
the quantum potential can be calculated explicitly [SI [23] but these complicated 
expressions do not give much to our intuition and we omit them here. (One 
can attempt at making the calculations by means of the tangent paraboloid 
technique [5] which was so useful above, but in the general case more powerful 
geometric methods [8j [23] are more safe and easier to implement.) Note that 
the Prokhorov method [30] can give the quantum potential too [2] but it is 
completely insensitive to the gauge structures. And let us finally mention that 
gauge structures (different from above) appear also in the algebraic approaches 
to quantization on a coset space [35]; we do not discuss it here. 

5 Conclusions 

We presented the main approaches to quantization of systems with second class 
constraints which do not involve path integrals: redefinition of the Poisson struc- 
ture (Dirac brackets), conversion of the constraints into the first class ones by 
introduction of new degrees of freedom and the thin tube approximation identi- 
cal to a way of realizing the holonomic constraints in classical mechanics. These 
methods give different results but often they involve very similar geometric con- 
structions and conditions. Many of the aspects of theory deserve a better under- 
standing. Geometric properties of the Abelian conversion are still unclear. And 
it's worth trying to describe the main features of higher codimensional Dirac 
quantization from an explicit operatorial perspective. This investigation can 
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be very important for our understanding of constrained quantum mechanics in 
general. 

A special remark is in order concerning the reference list. We do not attempt 
at making it absolutely complete. But, together with the references in the 
references, it should suffice to give a more-or-less full picture of research in the 
subject, at least as much as we were able to find it out in the literature. We 
think it is quite important because many similar results are scattered in different 
works being completely disconnected from each other. Jensen and Koppe [B] 
didn't know about the article of Marcus 29J, and neither of these works has 
influenced the results of da Costa [7J- The second article by da Costa [33] is 
rarely cited (as well as the pioneering one due to Marcus )• The authors of 
[10] and [11] were unaware of earlier works [5] and [13j . And in our articles [TJ[2] 
we didn't refer to very important works [5[ HI EH [28l [29l [33] due to the lack of 
knowledge. 
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